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I. INTRODUCTION 



There is growing interest in quasielastic electron scattering experiments with momen- 
tum transfers of order one GeV/c or larger. Clearly these require a relativistic treatment. 
Indeed, there have been many relativistic mean field and RPA calculations of quasielastic 
scattering These calculations often have large scalar S and vector V mean fields 

which shift the mass of a nucleon from M to M* , 

M* = M + S . (1) 

This changes the position of the quasielastic peak from u ^ q^/2M (with q the momentum 
and uj the energy transfer) to 

UJ ^ y^q2 + (M*)2 - M* . (2) 

This change agrees well with the "binding-energy" shift seen in a variety of experiments at 
moderate momentum transfers |q| < 600 MeV/c. 

However, Eq. (H) predicts very large shifts for momentum transfers of order 1000 MeV/c 
or larger. Such large shifts are not seen in data One would like to understand this 
limitation and develop a simple relativistic mean field model that is applicable over a broad 
range of momentum transfers. Indeed, most of the relativistic models assume the self- 
energies or mean fields {S, V) are independent of energy or momentum. This approximation 
is probably good at low momentum transfers but fails as the momentum transfer increases. 

Dirac optical model fits to proton-nucleus elastic scattering produce scalar and vector 
optical potentials which have much less energy dependence than the nonrelativistic mean 
field. When the Dirac equation is reduced to a Schrodinger like form the effective nonrela- 
tivistic mean field is 

Here E is the total energy of the nucleon (Tiab + M). The linear energy dependence for Uopt 
implied by Eq. (|]) agrees well with data for Tiab of order 200 to 300 MeV or less. However, 
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Eq. d^) (with constant 5", V) greatly overpredicts the energy dependence of Uopt above 300 
MeV. As a result Dirac optical fits need energy dependent S and V in order to reproduce 
elastic data at higher energies. One such fit is shown in Fig. |1|. 

The magnitudes of S and V decrease with energy. This trend is also reproduced in 
relativistic Brueckner calculations [0,3. In this paper we use such momentum or energy 
dependent potentials to calculate inclusive quasielastic electron scattering over a large range 
of momentum transfers. Results for a simplified form of momentum dependent potentials 
were presented in an earlier paper Our formalism is described in Sec. ||. Results for 



^^Fe and "^^Ca are given in Sec. IITII This section also discusses current conservation. We 



summarize in Sec. IV. 



II. FORMALISM 

The inclusive (e, e') cross section for excitation energy u and three momentum transfer 
q is 

2„ ,n4 /n2 



0"M 



^Rl{^,uj) + + tan2^)i?^(q,^)] , (4) 



dndE 

with the Mott cross-section cta/ = a^cos'^{9/2)/4:E'^sm^{9/2). Here = -q^ = - cj^ 
and 6 is the scattering angle. The longitudinal Rl and transverse Rt response functions are 
calculated in a local density approximation from nuclear matter at a density p = 2/c|,/37r^, 

i?^ = -— im(zngo + A^n^o) , (5) 

vrp 

RT = -±lm{ZU^,, + NU^^,) , (6) 
np 

for a target with Z protons and neutrons. (Note, q is assumed to be along the i axis so 
the subscript 22 refers to a transverse direction.) 

The polarization 11 is calculated with the nucleon Greens function G{p) and the electro- 
magnetic vertex F, 

n;.(g,o;) = -^| -^JAGiP + qKG{p)Tl] , (7) 
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for i = p (proton) or n (neutron). We assume the electromagnetic vertex has the simple 
form, 

Tl = Fl^, + Fi^ (8) 

even for off-shell nucleons in the medium. 

We calculate G in a mean field approximation. The nucleon self energy is assumed to be 

E = 7oF(p) + Sip) , (9) 

where 5* and V are taken from a Dirac optical fit at the self-consistent energy Ep, 



Ep = ^p' + {M + S{pW + V{p) (10) 
The Greens function for a Fermi momentum kp is 

r 1 77T 1 

Gip) = ip;Y + M;) ^ . + —Sipo - Ep)Q{kp - IpI) . (11) 

Here, 

Pl = {Po-V{p)M , (12) 

m; = m + s{jp), (13) 



and El = ^p2 + . 

We use S{p) and V^(p) from the real parts of an optical fit to p-'^'^Ca elastic scattering. 
This fit reproduces experimental data from 21 MeV up to a laboratory kinetic energy of 
1.04 GeV ip. For energies below 21 MeV we interpolate between this fit and the mean field 



approximation to the Walecka model values of |10 



Vmft = 354 MeV , (14) 
Sj^PT = -431 MeV . (15) 

For energies greater than 1.04 GeV we use the smooth extrapolation shown in Fig. |I]. (How- 
ever, results we present are not sensitive to this region.) Note, the optical fit has imaginary 



parts for S and V which we ignore. It remains to investigate the effects of these complex 
potentials. 

We use the fit potentials at the center of ^°Ca. However we scale these results by a, 

« = = (230/257)3 , (16) 

to represent an average over the nuclear volume. Here p(r) is taken from a relativistic mean 
field calculation of ^°Ca [jlO|. We assume the potentials scale with the density and then use 



a Fermi momentum of 

kp = 230 MeV/c , (17) 

in our calculations rather than a typical nuclear matter value kp ^ 257 MeV/c with po = 
2kp^ /Stt'^. This represents an average over the central and low density surface regions. 

The polarizations in Eq. (|^) are calculated numerically using the Greens function in 
Eq. ( pH]) as discussed in the appendix. The parameters of the calculation are the Fermi 
momentum kp, a and the two functions 5'(p) and V{p) from the Dirac Optical fit. This 
yields a simple relativistic impulse approximation that will be applied to data over a large 
range of momentum transfers in the next section. 



III. RESULTS 

In this section we present results for three calculations. The first is a free Fermi gas 
with ^(p) = V{p) = 0. The second is the Mean Field theory (MET) with momentum 
independent V and S given by Eqs. ([T^), (|15D (scaled by the a of Eq. (|16D). Finally, we 
consider momentum dependent ^(p) and V^(p). 

The longitudinal response for ^^Fe at |q| = 550 MeV/c is shown in Fig. ^ The free Fermi 
gas somewhat underestimates the average excitation energy of the response. This is better 
described by either the full momentum dependent calculation or the MFT. Indeed at this 
relatively low excitation energy of about 200 MeV there are not too large differences between 
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the momentum dependent S{p), V{p) and the MFT potentials. Therefore, either the MFT 
or the momentum dependent calculation can describe the average excitation energy of low 
momentum transfer data. 

We have used a simple local density approximation. This is expected to reproduce quali- 
tatively the position of the quasielastic peak. Full relativistic finite nucleus calculations (for 
momentum independent potentials) ref. agree well with the local density peak positions 
and do a good job of reproducing the detailed shape of the response. Note, the area under 
the longitudinal response is controversial. It is subject to systematic experimental errors 
and Coulomb distortion corrections. However, there is general agreement that (near this 
momentum transfer) there is a substantial binding energy shift compared to a free Fermi 
gas. This is seen either in separated longitudinal and transverse responses or in unseparated 
cross sections. 

Next, Fig. ^ shows the longitudinal response at a momentum transfer of |q| = 1.14 
GeV/c. At this momentum transfer the MFT response is at substantially too high an exci- 
tation energy. In contrast, either the momentum dependent or the free Fermi gas calculation 
provides a reasonable description of the position of the quasielastic peak. Thus the reduction 
of 5'(p) and V{p) with increasing energy or momentum (as shown in Fig. |I]) corrects (at 
least in large measure) the tendency for the MFT to overpredict the binding energy shift. 
At still higher momentum transfers, we expect small differences between the momentum 
dependent calculation and a free Fermi gas because S{p) and V{p) continue to decrease. 

The transverse response for ^''Ca is shown in Fig. ^ for |q| =550 MeV/c. Again, at this 
q the MFT and momentum dependent calculations are similar. Note, there is substantial 
extra strength in the transverse response at high excitation energies from meson exchange 



currents and Delta production, see for example . These are not contained in our simple 
model. 

Note, a free Fermi gas reproduces the peak position well at high q. However, there is no 
theoretical motivation for simply ignoring the mean fields which are known to be present. 
Furthermore, a free Fermi gas will fail at low q where it is generally believed one needs a 



binding energy shift. Thus, we think, the only theoretically consistent model is one which 
includes mean fields and takes into account their momentum dependence. 

Alternatively, some nonrelativistic models simply assume a momentum independent exci- 
tation energy shift equal to the average binding energy of a nucleon. This works reasonably 
well at low q and possibly could be extended to high q by using some form of relativis- 
tic kinematics. However, there is no theoretical justification for this prescription. Strictly 
speaking one would have to assume an energy dependent potential which is big for bound 
states and then goes rapidly to zero for positive energies. This is unrealistic. Thus, we 
believe, a theoretically consistent description will have to address explicitly the momentum 
or energy dependence of the mean fields. 

We now discuss current conservation. The full electromagnetic vertex should satisfy the 
Ward-Takahashi identity, 

q,T%,, = G-\p + q) - G-\p) . (18) 

This is no longer true for the vertex F^, of Eq. given that G is calculated with momentum 
dependent potentials (and that Eq. (|^) has form factors). Clearly we must add a vertex 
correction. However, there is no unique way to calculate this correction without detailed 
knowledge of the processes included in S'(p) and V{j>) and in Fi and F2. 

Instead, we add a very minimal vertex correction so that not Eq. (0) but at least current 
conservation is satisfied, 

[/(^(p), V(p))g^r^„,[/(^(p), y(p)) = . (19) 

Here, the spinors U satisfy a Dirac equation with S", V and 

r^,,, = r'^ + Ar^ . (20) 

We choose, 

^T^' = q'^^iSip + q) + loVip + q) - ^(p) - 7oV^(p)] • (21) 
Qx 
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This simple prescription insures, 

q^U^, = U.^q'' = . (22) 

Calculations using Eqs. (0), (l2l|) are shown in Figs. H, ^and are only slightly smaller than 
calculations with AF^ = 0. Therefore, we do not expect large errors from our (slight) 
violation of current conservation. However, Eq. (^) is non-unique and this point deserves 
further study, see below. 

IV. SUMMARY AND CONCLUSIONS 

There have been many relativistic mean field calculations for inclusive quasielastic elec- 
tron scattering. Ironically, although these do a reasonable job at low to moderate momentum 
transfers, many of these relativistic calculations are not appropriate for momentum trans- 
fers near one GeV/c or above. This is because they assume the mean fields are independent 
of momentum or energy. As a result they predict too large a binding energy shift in the 
position of the quasielastic peak. 

Instead, we assume a simple model where the momentum dependence of the scalar and 
vector mean fields are taken from an optical model fit to p-^°Ca elastic scattering. The 
inclusive response is calculated in a relativistic impulse approximation. Reasonable predic- 
tions are made for the position of the quasielastic peak at both low and high q. Therefore, 
the model can be used over a large range of momentum transfers. 

Our calculations are for nuclear matter in a local density approximation. We expect 
this to qualitatively reproduce the position of the quasielastic peak but not its detailed 
shape. One should repeat the calculations in a full finite nucleus formalism, such as from 
reference 0, using energy dependent potentials. 



^Note, Eq. ( ^l[ ) does not contribute to Rl when dotted into the conserved electron current. 
However, Eqs. (^), (^ assumed current conservation. One gets a different answer for Rl from 
Eq. (^) if one sets AF^ = because in this case the current is not conserved. 
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We have only briefly mentioned vertex corrections. These are needed for current con- 
servation in momentum dependent mean flelds. Future work should calculate the vertex 
correction in a full microscopic model. For example, one can choose a (somewhat unusual) 
set of meson couplings (including nonlinear self-couplings) so that a relativistic Hartree-Fock 
(HF) calculation approximately reproduces the momentum dependence of the mean flelds. 
This allows the vertex correction to be calculated explicitly from the diagram in Fig. ^. 
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APPENDIX 



Here we summarize the calculation of imaginary part of the polarization 11^'^. The Greens 



function of Eq. (11) can be written in three parts; (1) the propagation of an antinucleon in 
the Dirac sea, the propagation of a (2) hole inside and a (3) particle outside of the Fermi sea. 
Since vacuum polarization does not contribute to the imaginary part for the space-like g^, 
the propagation of antinucleons can be dropped and the Greens function Eq. (0) written, 



Gip) 



2E; 



e{kF-\p\)_^ e(|p|-M 



(Al) 



_Po - Ep-ie po- Ep + ie_ 

Using this Greens function, it is straight forward to integrate over po and the imaginary part 
of the polarization becomes, 

p^(i|p| dcos6 



im n 



liu_ 



X [Tr + T^" + Tg^l 



-e(|p + q| - kp) Qikp - IpI) 5iEp+^ -Ep- go) 



(A2) 



where cos6' is the angle between p and q. Here Ti,T2 and T3 are from the trace of Dirac 
matrices, 



^ 2M 

,2 r 



rr = [{p + qrp*" -p*-{p + gyg^-" + p*^{p + qr + M;M;^^g^^] , 
m; [{p + qr-q g^" -{p + g)*^g1 + m;^^ [p*y^ - r ■ qg 

[M;^^M; + p*-{p + q)*] {g^^ ql - q'^q^ 
+p* ■q[{p + q^q'' + {p + qTq^] - ql [{p + g)*V^ + {p + qVp*'. 
Hp + qy-q W + -2p*-q{p + q)* ■ q g'^-' 



rp/liy 
-L ^ 



4M2 



(A3) 
(A4) 



(A5) 



with p*^ = {Ep — V{p), p}. Equation ( |A2| ) is evaluated numerically. The phenomenological 
self-energies (l^(p), ^(p)) are fitted to polynomials in p. Finally, the angle integration of 
the (5-function restricts the momentum integration so that |cos0| < 1. The form factors F^ 



and F2 are taken from Ref. in . 
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FIGURES 

FIG. 1. Scalar S and vector V mean fields versus laboratory kinetic energy Tiab- Solid lines are 
real parts of a Dirac optical model fit to p-^'^Ca elastic scattering ^ while dot-dashed lines are 
the constant relativistic mean field values [|l^ . The dashed lines are extrapolations as described in 
the text. 

FIG. 2. Longitudinal response function for ^^Fe at a momentum transfer of q=550 MeV versus 
excitation energy tu. The solid curve is the full momentum dependent calculation including the 
vertex correction from Eq. (|2l| ) while the dot-dashed curve omits this vertex correction. The 
short-dashed curve is the response for a free Fermi gas while the dashed curve is the response 
assuming constant relativistic mean field theory self-energies. The data is from Ref. 

FIG. 3. Longitudinal response for ^^Fe at a momentum transfer of q=1.14 GeV/c. See the 
caption to Fig. The data is from Ref. @. 

FIG. 4. Transverse response for ^''Ca at a momentum transfer of q=550 MeV/c. The solid 
curve is the full momentum dependent calculation while the short-dashed is a free Fermi gas and 
the dashed is for constant relativistic mean field theory self energies. The data is from Ref. 

FIG. 5. Vertex correction in a relativistic Hartree-Fock approximation. The nucleon line is 
solid while the dashed line is an exchanged meson and the wavy line a photon. 
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